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Abstract 

Coulomb corrections for quasi-elastic scattering of electrons by nuclei are calculated using eikonal distorted 
waves. Corrections to the lowest-order eikonal approximation are included in order to obtain accurate results. 
Spin-dependent eikonal phase shifts are evaluated and they yield very small corrections to the longitudinal 
and transverse cross sections at electron energies of 500 MeV or higher. Because of this the Rosenbluth 



procedure is accurate for separation of the longitudinal and transverse response functions. The effective- 
momentum approximation is also found to be accurate with regard to removal of the remaining Coulomb 
effects from the distorted waves. Calculations are presented for electron scattering from 208 Pb and 56 Fe 



■ momentum approximation is also found to be accurate with regard to removal of the remaining Coulomb 

nuclei at energies of 500 MeV and 800 MeV and momentum transfers q — 550 MeV/c and 900 MeV/c. 
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I. INTRODUCTION 



Measurements of cross sections for quasi-elastic electron scattering provide basic tests of our 
understanding of nuclei. Experiments have been performed at the MIT Bates Laboratory ] , [], 
, at the Saclay Laboratory [k], [n], [12} Q, Q, Q] and at SLAC Q, Q, H] in 
order to explore this reaction. A recent review of quasi-elastic scattering provides a guide to the 
experimental and theoretical results [19]. Generally it is assumed that the reaction is dominated 
by the exchange of one hard photon between the electron and a nucleon in the nucleus. Many 
soft photons also are exchanged and their effects are described by distorted waves based on the 
Coulomb potential in the initial and final states of the electron. An important issue is to account 
for the effects of the Coulomb interaction in a manner that allows extraction of the nuclear response 
functions, Rl and Rt, which correspond to plane- wave matrix elements of the longitudinal and 
transverse parts of the currents. After the nucleon form factor is divided out, the longitudinal 
response function at a fixed and sufficiently large value of the momentum transfer, q, is expected 
to satisfy a Coulomb sum rule, i.e., J duSi^q,^) ~ Z, where Z is the number of protons. Such 
a sum rule should hold under general conditions for a nonrelativistic description of nuclear wave 
functions and currents provided only that q is sufficiently large to make the effects of correlations 
small. Although the response functions are affected by the final-state interactions (FSI), the sum 
rule should not be affected by them because it relies on a sum over a complete set of final states 
of the nucleus. When relativistic effects in the nuclear current are taken into account, there are 
minor relativistic corrections to the sum rule but it remains useful. However, attempts to verify it 
experimentally have produced puzzling results and controversy about Coulomb corrections. 

In order to obtain the longitudinal response function, it is conventional to perform a Rosenbluth 
separation after dividing the cross section by the Mott cross section om and some kinematical 
factors, 

with Q 2 = uj 2 — q 2 and where om = 4:a 2 E 2 cos 2 (^6 e )/Q 4: . The longitudinal virtual-photon polar- 
ization is defined as 

At fixed q and uj, it varies from to 1 as the electron scattering angle 8 e varies from 180 to 
degrees. Measurements of the cross section at different values of e for fixed values of q and uj allow 



a separation of the experimental response functions, R e ^ pt and R^ pt . The R e ^ pt that is extracted 
from experimental cross section as in Eq. ([I]) can differ from the desired Rl, which is the plane- 
wave matrix element of the longitudinal current, for two reasons. One is that Coulomb corrections 
associated with spin-dependent effects in the electron wave functions can cause the contributions of 
the longitudinal and transverse currents to have different dependence on 6 e than is assumed in the 
Rosenbluth separation. [^(J The other is that Coulomb distorted waves alter the response. These 
Coulomb corrections should be removed before the Coulomb sum rule can be evaluated. In addition, 
nucleon form factors should be divided out of the longitudinal response function. A standard 
approach to modeling the Coulomb corrections is to use the distorted wave Born approximation 



(DWBA). One solves the Dirac equation for the electron distorted waves in t 
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re presence of the 

" n 



22j, 123, l2J When 



Coulomb potential for both the initial and final states of the electron, 
these distorted waves are combined with nuclear wave functions obtained from a model of nuclear 
structure, and appropriate current operators, cross sections may be calculated and compared with 
experimental results. The DWBA analysis involves extensive numerical calculations. Owing to the 
infinite range of the Coulomb potential, partial-wave expansions converge very slowly, the more so 
as the energy increases. 

A number of works have used the eikonal approximation in order to simplify the analysis at high 
energies. The use of the eikonal approximation also must be accompanied by inclusion of "focusing 
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27( Czyz and Gottfried [281] used 



factors" that are not present in the eikonal wave functions, 
the eikonal approximation to analyze electron scattering but that work did not include focusing 
factors. An analysis based on the Schrodinger equation showed that the corrections to the eikonal 
approximation generally produce a focusing factor in the wave function. 
al. also is based on the eikonal approximation 
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29( Work by Giusti et 



311 ] and some recents works have combined the 



eikonal approximation with semi-classical focusing factors in order to assess Coulomb corrections 
in quasi-elastic scattering. 
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33] 



A very simple effective-momentum approximation (EMA) for treating the Coulomb corrections 



was developed by Rosenf elder 



27( and Triani et al. [3J, |35j. In the EMA the effects of the Coulomb 



potential are incorporated as shifts of the initial and final electron momentum values that should 
be used in a plane-wave Born approximation (PWIA) analysis. The shifted electron momenta are 
the effective momenta. They imply a corresponding shift of the photon momentum, q — * q e ff- 
Interpretations of quasi-elastic data depend upon many experimental details and different ex- 



periments have produced significantly different values of the Coulomb sum rule 



. HQ 



In addition 



to possible experimental differences, there are theoretical differences in the analysis of the Coulomb 
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corrections because the DWBA analysis based on partial waves has been used for some experiments 
and the EMA analysis has been used for others. Sometimes it is assumed that the nucleon form 
factors can be pulled out of the matrix element and evaluated at the momentum transfer of the 
electron. [3] At other times the form factors are evaluated at the effective photon momentum, 

In order to address questions about the theoretical differences in the treatment of Coulomb 
distorted waves, we developed a systematic eikonal expansion in Ref. [37] that provides more 
accurate eikonal wave functions for a DWBA analysis. The accuracy is good enough to eliminate 
concerns about use of the eikonal approximation at the energies of interest. Moreover the "focusing 
factors" arise naturally as part of the corrections to the eikonal approximation and the ad-hoc 
procedure of incorporating them is replaced by a systematic procedure. The eikonal expansion was 
found to converge rapidly at electron energies of interest. It has the advantage of providing insight 
into the nature of the Coulomb corrections because the focusing factors, the eikonal phase shifts 
that determine the momentum shifts and spin-dependent effects can be isolated for study. 

In order to assess the accuracy of the EMA using the eikonal wave functions, we used in Ref. 37] 



a very simple model of the nuclear wave functions and we neglected the spin-dependent Coulomb 
corrections for simplicity. In order to compare full DWBA calculations of Rl with the results based 
on the EMA, the DWBA results were fit to the EMA formula, 

R L (q,u;)=AR^ WIA ( qeff ,u J ), (3) 

where the effective photon momentum is given by 

q e jf = ki ki — 5k — kj kj — 5k . (4) 

There are two parameters in our EMA fits: the momentum-shift 5k and the overall normalization 
constant A ~ 1. Usually the parameter A is assumed to be unity when experimental data are fit 
using Eq. ©. In model calculations that assumption can be checked because the PWIA response 
is known in the model. A factor /ema can be used to relate 5k to the Coulomb potential at r = 0, 
as follows, 

fEMA = W)- (5) 

The factor /ema * s approximately the same for different nuclei. We found momentum shifts in 



Fig. 7 of Ref. 



371 ] that correspond to /ema(^) ~ 0.7 near the peak of the response function for 



500 MeV e scatttering. Larger values of Jema up to about 1 were found at the smallest and 



largest to values but the response is small at those points. Although fits of the DWBA results 
can be made more precise by allowing 5k to depend on the energy loss, u, in this work we use a 
constant shift 5k. That yields reasonable results and is simpler and thus preferable for the analysis 
of experimental data. 

If the momentum-shift 5k and normalization constant A are determined theoretically for a given 
nucleus such that the DWBA response is well described by the EMA fit of Eq. ([3|), then one may 
equate i?i(q, lu) at fixed electron beam energy, E, and fixed momentum transfer q to a constant 
A times the PWIA response function evaluated at the effective momentum transfer. That would 
remove the Coulomb effects to a reasonable approximation and allow the PWIA response to be 
extracted from experimental data. We expect similar Coulomb corrections for a variety of nuclear 
models. [39] The goal is to remove them with minimal reliance on any nuclear model. However, it 
must first be determined how well Rl can be extracted from experimental data. 

The spin-dependence of the eikonal wave functions was omitted in our previous paper, which 
left unanswered the question of the accuracy with which the desired response functions might be 
extracted from experimental cross sections. That is the first issue addressed in this paper. In 
Sec. II we restate the essential results of the eikonal expansion for Dirac wave functions and focus 
on the spin-dependent eikonal corrections. These are shown to provide very small differences to 
quasi-elastic cross sections, i.e., the helicity matrix elements of the electron current are very close 
to those based on the PWIA. The consequence is that when both initial and final electron energies 
are 200 MeV or more, the usual Rosenbluth separation provides an accurate separation of Rl and 
Rt, well within the limits of experimental accuracy. Although the Rosenbluth separation should be 
accurate, there remain significant Coulomb effects within Rl- They can be treated with reasonable 
accuracy by use of the effective momentum approximation. 

In order to determine more realistic values of the momentum shift, 5k, shell model wave func- 
tions are used to describe the nucleus in this work. Our calculations are simplified by use of 
an approximation that is introduced by us in Ref. 38|] and described in Section III, and which 



is denoted EMAr. That approximation applies the effective momentum approximation to the 
hard-photon propagator and form factors in order to reduce the numerical evaluation to a three- 
dimensional integration that provides a careful treatment of the full r-dependence of the Coulombic 
effects from the electron wave functions. In the EMAr analysis, the nuclear current is handled in 
terms of a hadronic tensor, which can be extended to include the neutron contributions to cross 
sections. A comparison of the full DWBA and EMAr calculations for the Is shell of 208 Pb shows 
close agreement of the results. 
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Section IV presents numerical calculations for quasi-elastic scattering by 208 Pb and 56 Fe using 
kinematics that are relevant to a recent experiment at Jefferson Laboratory. The nuclear model used 
is simple and a number of significant effects are omitted from the calculations, such as final-state 

but the calculated cross 



interactions 



40l | . correlations, 



41 



43 



42( | and pion and A production 
sections are expected to be roughly similar to experimental ones. The main goal is to determine 
suitable fitting parameters for use in applying Eq. ([3j) to experimental data so as to determine 
R^ WI . Conclusions are presented in Section V. 



II. QUASI-ELASTIC RESPONSE FUNCTIONS 

Because electron energies of interest generally are much greater than the electron mass, and the 
Coulomb potential and photon exchange are vector interactions, electron helicity is conserved to 
a very high degree of accuracy in quasi-elastic scattering. For example, using 500 MeV electrons 
one finds that the helicity is conserved except for terms of relative order m 2 /E 2 ~ 10~ 6 . In this 
work we keep only the effects that are consistent with helicity conservation. 

The distorted-wave Born approximation is used with eikonal wave functions for the electron. 
For outgoing-wave (+) or incoming-wave (-) boundary conditions, the Dirac wave functions for 
potential V(r) are written as [371] 

\2Xu^(r)J 

«C±)(r) = fl-^Y^V^V^V^, (6) 

where £\ is a two-component helicity spinor and A = ±| is the helicity eigenvalue. The lower com- 
ponents of the Dirac spinor are simply 2A times the upper components because the electron mass 
is neglected. The wave propagates in the z-direction, which is along the asymptotic momentum k, 
and an impact vector b is defined as the part of r that is perpendicular to the 5-direction. The 
eikonal phases % ■> ^ an d 7^ = 7^ ± £<$w are obtained from integrals over the potential 



along the z-direction as shown in Ref. 
the energy is E and E% = E + m. 



371 ] . The spin matrix in the eikonal phase is a e = a ■ b x z, 
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A. DWBA analysis 



The DWBA cross section for knockout of a nucleon of momentum p involves a two-dimensional 

integration over the angles of the knocked-out nucleon as follows, 

1 da Q 4 ^ f pE p | ,2 

x ' dn p —^\M n i m \ ■ (7) 



a M dndw C0S 2^)^J p (2vr)5 

Omitting the final-state interactions of the nucleon, the matrix element for quasi-elastic knockout 
involves a six-dimensional integration, 

= / ^ J Ae'^'h^rlJ/^ ^^^V n^tq' " P) (8) 

where jA^(q', p) is the nucleon current, he is a four- vector of helicity matrix elements of the electron 
current and ipnimi 1 ^ — p) is the momentum-space wave function of a nucleon in the nucleus with 
quantum numbers n (radial), 1 (angular momentum) and m (z-component of angular momentum). 
Subscripts i and / refer to the initial and final electron states with asymptotic momenta kj and kj 
that provide the respective z-directions for incoming and outgoing waves. The exchanged photon 
has energy u) and momentum q', the initial electron helicity is Aj and the final electron helicity is 
Xf. The momentum transferred by the electron is q = kj — k^ and it differs from the momentum 
q' of the photon because of the Coulomb effects. Nucleon form factors within the nucleon current 
depend on the photon momentum, q', that is integrated. The sum of eikonal phases for incoming 
and outgoing waves is \ = X) ^ ( r ) + Xi + ^ ( r ) an d the focusing factors are defined by 

fi{r) = (l ~ 

//W = (l- — ) (9) 

where E<n = Ei + m and E^f = Ef + m. 

The Coulomb potential for scattering from a 208 Pb nucleus is shown in Fig. ll The solid line 
shows the potential based on a fit of experimental data for the charge density 4a] and the dash 
line shows the simple potential used in our calculations, namely, 

".M = l^TWn- < 10 > 

The simple potential allows an analytical evaluation of the eikonal phases. The parameters Vq 
and R are chosen so that the simple potential has the same value as the empirical potential at 




r = and the same average value as the potential based on experimental data in the sense that 
f dr p expt (r)V c (r) = J dr p expt (r)V exp t(r) . This ensures that V c (r) provides a good fit in the range 
where the nuclear density is significant (the dot-dash line shows one-tenth the nuclear charge 
density for reference). Eikonal phases based on the Coulomb potential of Eq. (fTU|) are shown in 



0.03 




r (fm) 



FIG. 1: Coulomb potential for Pb. The solid line shows the potential based on the empirically determined 
charge density. [45J The dash line shows the approximate potential used in this work based on parameter 
values Vo — 0.0256 GeV and R —7.1 fm. The dash-dot line shows one-tenth the charge density, 0.1p c / 1 (r), 
for 208 Pb. 

Fig. [2 The phases x[ + ^ an d oj^ ^ are shown as a function of impact parameter |b| at z = and 
the nucleus is located at b = 0, z = 0. Note that Xi + ^ anc i X) ^ each have the same behavior at 
z = and the total phase x = Xi + X/ 18 approximately double the values shown. Figure [2] 
shows that the eikonal expansion produces well converged results for the wave function at electron 
energy equal to 500 MeV for scattering from 208 Pb. 
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FIG. 2: Eikonal phases at z = for 500 MeV electron scattering from m Pb. 

Coulomb corrections affect the helicity matrix elements of the electron current because of the 
spin-dependent eikonal phases involving as follows, 

Kiv) = 5 A/Ai Ci / (^)e i<Je ^/" ), '{l ! 2A i a} e ^# ) eA i . (11) 

I — I 

As shown in [37|, the required helicity matrix elements are given by 

h° e = A 2Xl cos^9 e + C 2 \ i sin-9 e — > cos^9 e , 

K = B 2Xi sin^9 e + D 2Xi cos^9 e — > sin-9 e , 

h\ = 2i\i(^A 2 \ t sin^9 e - C 2 x t cos^9 e ^ — ► 2iXisin-9 e , 

K = B 2Xl cos^9 e - D 2Xi sin^9 e — > cos^9 e . 

(12) 

The Rosenbluth separation implicitly assumes that the helicity matrix elements take the plane- wave 
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values that are shown following the arrows in Eq. (|12p . One sees in Eq. (|12p that the spin-dependent 
Coulomb effects cause each component of the four-vector of electron helicity matrix elements to 
involve both cos^9 e and sin^9 e . This property carries over to the longitudinal and transverse 
currents. Consequently, R e ^ pt and R^ pt extracted using the Rosenbluth separation each involve 
admixtures of the longitudinal and transverse currents. 

The spin-orbit parts of the eikonal phases enter the helicity matrix elements in the following 
four combinations as shown in Ref. 

A 2Xi = cos^f ] *cos^ +) - sin^sm^'V 2 ^*^', 
B 2Xl = cos^f ] * 'cos^ +) + sin^sinj^e 21 ^-^, 



C 2Xi = cosrff ' sinj> ' e ZlAi<pl + cos^y 1 siwyK 
D 2Xi = cos^f ] * sin^\ +) e 2tX ^ - cos^ +) ainrff^ e~ 2iAi ^ . 

(13) 

Figure [3] shows the eikonal phases 7^ + ^ and 5^ that are the real and imaginary parts of 7W . 
There is a simple relation between the phases as follows, 

- (±) (r) = 7 (±) (r) ± l6 (±) (r) = _}_ d_ ^ (±) (r) ± . J±) (r) ^ (14) 



This relation has been corrected from the one given in Ref. because a factor ^ was omitted 
there. 

Our previous numerical calculations omitted the spin-dependent Coulomb corrections, using 
instead the PWIA helicity matrix elements that are indicated following the arrows in Eq. (|12p . 
In this work we include the spin-dependent Coulomb corrections and find the important result 
that they provide negligible corrections to quasi-elastic cross sections. This can be understood 
qualitatively as follows. For a 500 MeV electron scattering from the Coulomb potential of a 
208 Pb nucleus, the phase 7 ( - + * ) = 7^ + i5^ typically has magnitude of 0.03 or less as shown in 
Fig. O It vanishes at zero impact parameter. Consequently, the coefficients A 2Xi ~ 1 + C(7) 2 and 
E>2\i ra 1 + C(7) 2 are very close to unity. The coefficients C 2Xi and D 2X . have contributions that are 
first order in 7. They also involve phase factors e 21 ^'^ or e 2lXi ^f , where tanfif = tanfii/ (^cos^ — 
sin^f / (tan0icos<f)i)^ with 6i and (f>i being the angles of the vector r relative to the polar axis along 
ki. When current matrix elements are integrated over r and summed over helicities, cancellations 
stemming from these phase factors cause the contributions to be very small. Current matrix 
elements are Fourier transforms and the faj dependent phase factors can receive support in the 
integration. However, the quasi-free cross sections involve an additional integration over angles of 
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FIG. 3: Spin phases for 500 MeV electron scattering from 208 Pb. 

the knocked-out proton. The net effect is to reduce substantially the contributions from terms that 
involve the azimuthal angles <j>ij- Our numerical results show that cross sections calculated with 
the spin-dependent Coulomb corrections included are closer than one part per thousand to ones 
calculated using the PWIA helicity matrix elements. An example of this is shown in Figured! which 
shows the ratio of longitudinal response functions (calculated from Eq. (|30p ) with and without the 
spin-dependent Coulomb corrections at 500 MeV electron energy. Only for large energy loss, where 
the final state electron energy becomes small, does the ratio differ from unity by more that a few 
parts in ten thousand. For reference, the dotted line shows the variation of 1 + .OOOIRl- Note 
that for energy loss u) = 300 MeV, the final electron energy is 200 MeV, but the ratio of cross 
sections with and without the Coulomb spin corrections differs from unity by less that one part 
per thousand. As the energy loss increases, the Coulomb spin corrections become relatively more 
important but the response function is decreasing to zero. The net effect is that the absolute error 
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in the response function is about one part in ten thousand of the maximum value of Rl- This 
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FIG. 4: Circles show the ratio of the longitudinal response function R^ that includes the spin-dependent 
Coulomb corrections to the longitudinal response function that omits them for 500 MeV electron 
scattering from 208 Pb at q=550 MeV/c. The dotted line shows 1 + .0001i? L . Note that the effects of the 
spin-dependent Coulomb corrections are largest for uj- values where the response function is a small fraction 
of its maximum value. 

is important because the helicity matrix elements govern the dependence of the cross sections on 
the electron scattering angle, 8 e , at hxed q. When the PWIA helicity matrix elements provide 
an accurate approximation, the angle dependence is the same as in PWIA cross sections, and 
the Rosenbluth separation can be used to extract the longitudinal and transverse current matrix 
elements. We find negligible mixing of actual longitudinal and transverse current matrix elements 
in the response functions extracted using the Rosenbluth separation for all the cases evaluated in 
this paper. 

Significant Coulombic effects remain in the Rl and Rx response functions that can be extracted 
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by use of the Rosenbluth separation because of the spin-independent Coulomb effects. They may 
be treated using an effective momentum approximation. 



III. EFFECTIVE PHOTON MOMENTUM APPROXIMATION 

Because the effective-momentum approximation was found to be accurate in our prior work, in 
this work we approximate the photon momentum q' that appears in the photon propagator and 
the nucleon current of Eq. ([8]) by an effective momentum q e // as follows, 

(g. q') jfip^ff) 

q' z -w 2 q^-w 2 

where q e / f is given in Eq. ([4]) . Note that this approximation also evaluates the nucleon form factors 
within the nucleon current at the effective photon momentum. This is a minimal use of the effective- 
momentum approximation designed to reduce the computation to a three-dimensional form, i.e., 
the approximation allows the photon propagator and the nuclear current to be factored out of the 
integral over q', which is then performed to obtain J e~ t<l '' r tp n i m (q L ' — p) = e~ ip ' r ip n i m (r). This 
yields 

M nlm — M™ Ar = £-$(P, Qeff) M» lm (16) 
where Qgft = leff ~ ^ ano - we define 



^nim = J MQ-vy'jxU /i(r)//(r) K(r) ^ m (r). (17) 



This procedure has been discussed in Ref 38] and is called the EMAr approximation. It has the 
advantage over the usual EMA of not approximating the r-dependence of the eikonal phases and 
focusing factors that provide the Coulomb corrections to the electron wave functions. The three- 
dimensional integration of Eq. (|17p provides a good correspondence with the full DWBA analysis 
at much lower computational cost. Procedures to determine appropriate values of 5k, and thus the 
effective photon momentum that is factored out of the integral, are discussed further on. 

A. Hadronic tensor 

In the EMAr analysis, the bound-state nucleon's wave function is taken to be a product of a 
Dirac spinor, u(p — q), times a nonrelativistic wave function for a nucleon, i.e., (M/(£ , p _ q ) 1 / 2 u(p — 
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q)Vvum(p — q), and the knocked-out nucleoli's wave function is a Dirac spinor, (M / ' E p ) 1 I 2 u(y>). 
The relevant nucleon current is 



i N 



K^ 2 u(p) [ 7 „Fi + 2M^2^^J u(p - q), 



(18) 



where Fi(Q 2 ) and F 2 (Q 2 ) are nucleon form factors, k is the anomalous magnetic moment and 
K = M 2 / (EpEp-q) is a normalization factor arising from the spinors. The EMAr cross section for 
knock-out of a nucleon of momentum p by absorption of a photon then takes the form 

thiL - Jw) £ / dn 'w?w^ M " M w)m:l (19) 

where the hadronic tensor is, 

tf + M ( 



a \j>-i + M 



Carrying out the trace over nucleon spins produces 



(20) 



W^(p,q) = -g^F 2 + 



P^{p ~ l)v + (P~ Q)iiPv -p-(p- q)g fll y 2 



2M 2 



+ 



t?1 



,2M 2 ' 



- 



~8M 2 2 
-q 2 {2p ll p u - p 2 g ilv ) 



QiiQvip +P- q) +^{%Pv +P^qu){q -P- q) + {2p ■ q- q )p ■ qg^ 



(21) 



The hadronic tensor is gauge invariant when the momenta are on mass shell, i.e., p 2 = M 2 and 
(p ~ q) 2 = M 2 . These conditions require that p ■ q = \q 2 . We also use E p = M + uj for an 
initial nucleon at rest, leading to if = M 2 /[M(M + uj)]. Using the on-mass-shell kinematics and 
Q 2 = —q 2 leads to the gauge invariant form that is used in this work, 



W^(p, q) = —(p lx - -q^ (p v - -q u ^ (if 



k 2 Q 2 



+ 



Fi + 



^ixv - q^qu 



Fi + kF 2 



AM 2 " 2 ) ' AM 2 

This form of the hadronic tensor is evaluated at q^ — > (w,q e jj) in the EMAr analysis. 



\ 2 

) • (22) 



B. Cross sections and response function 



The form of the hadronic tensor shows that the cross sections involve an incoherent sum of two 
parts, which is a consequence of averaging over nucleon spins. The cross sections may be written 
concisely in terms of the Sachs form factors, 



G M = F 1 + kF 2 



(23) 
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using the combination of form factors, 



r 2 _i_ Q r<2 

E ANP^M 



1 + 



(24) 



in place of F- 2 + F$. Using current conservation, which takes the effective forms M. n im • 4e// = 
(uj/q e ff)Mni m and h e -<? e // = (u/q e f f)h®, the longitudinal components of the electron and nuclear 
currents can be expressed in terms of the correspond charge components. In so doing we arrive at 

2 



1 da Q 4 v f pE p K 

a M dn<ko cos\^)^J P (2vr)3(Q2 // )2 



M + \uj 



M 



-M 



nlm 



U) \ p - T 



Q2 



(25) 



The transverse amplitude arising from the vector part of the convection current is 



Mni m = M 



nlm 



Mnlm ■ q e // )q e // 



(26) 



As has been shown, it is a good approximation to omit the spin-dependent eikonal effects. Then 
the helicity four-vector simplifies to the plane-wave form, 



uuPW x f &e ■ &e r,\ ■ @e @e\ 

ivz =d\ r \ <cos — , sin — .ZAiSin — , cos — >. 



(27) 



It follows that 



\M 
\M 



nlml 



nlml 



(X COS 



oc sin 



nlm 



X 



K PW 9 n,K PW = -2W|. 



(28) 



The interference terms between M° n i m and -M^im van i sn by symmetry and the quasi-elastic cross 
section takes the form of Eq. (pTJ) with response function Rl being proportional to the square of the 
matrix element of the time-component of the current. The longitudinal cross section is obtained 
as 

2 



1 do L 



Q 



pE p K 

\3 (n 2 \2 



M+±LU 



M 



M 



o 

nlm 



1 



uJ 



<ff 



GUQiff)- (29) 



om dSlduj cos^f) ""'W (Ql ff ) 

This analysis suggests that to extract the longitudinal response function, one should perform 
the Rosenbluth separation as in Eq. ([!]) except that a factor ^ff/Qtff should be used in place 



15 



of the g 4 /Q 4 kinematical factor in order to cancel the 1 — LJ 2 /q^^ factor in the matrix element. 
Applying the factors cQft/Qttt anci 1/fAf to the longitudinal contribution of Eq. ([25]) yields 



u EMAr 



M(M + u) 



Q 4 



g\(q 2 



eff) 



2^J d ^ (27r) 3 

nlm 



M 



o 

nlm 



(30) 



The first prefactor arises from the relativistic wave function normalization factor, K, and the 
current, where p° = M + lu and q° = uj have been used. The second prefactor should cancel 
to a large extent with similar factors in the .M?, amplitude. In particular, the helicity matrix 
element h® ~ cos{^f) cancels to a high degree. To the extent that the focusing factors from 
the electron wave functions in the DWBA matrix element can be approximated at r = 0, they 
give approximately a factor 1 k . k f f — to the matrix element, which should cancel with Q 4 /Qfj^ = 

( kjkf 
\kf f k e / f 



from the prefactor. Because the EMAr analysis involves an integration over r with 
the full coordinate dependence, the focusing factors can differ from the approximate result. For 
example, the contributions from radial wave functions that vanish at r = 0, as in the shells with 
L > 0, depend on the focusing factors away from the origin. Consequently, accurate values of 8k 
are not known a priori such that the second prefactor cancels the Coulomb effects owing to focusing 
factors. A procedure is needed to determine them. 

In order to test the Coulomb sum rule, one should remove the nucleon form factor. Consistent 
with the photon momentum in the nuclear current being shifted, the form factors evaluated at 
Qeff should be divided out of the cross section, giving 

1 R L {q,u) 



z G%(Qi ff y 

dwSL{q,uj) 1. (31) 



C. Comparison of EMAr and DWBA analyses 

The six-dimensional (6D) integration of Eq. ([8]) should be performed for nuclear models in order 
to obtain the full DWBA response function. Similar calculations of the EMAr response function 
should then be normalized, by choice of the Q e ff value in the prefactor of the EMAr amplitude, such 
that the magnitude of the peak DWBA response is reproduced by the EMAr analysis. However, the 
full six-dimensional integration in the matrix element together with the two-dimensional integration 
over angles of the final momenta is extremely time consuming when many shells contribute to the 
response. In this work, a limited form of the full DWBA analysis has been performed based on 
the wave function of a single shell, the Is shell of the 20S Pb nucleus. We then normalize the EMAr 
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analysis to the DWBA analysis for the ls-shell. Figure [5] shows the PWIA, EMAr and DWBA 
response functions based on the Is shell. In the full DWBA calculation, the nucleon form factor 
F\q' 2 — a; 2 ), the photon propagator, l/(q' 2 — uj 2 ), and the current-conservation factor ^1 — u; 2 /q /2 ^ 
have been evaluated at the running photon momentum that is integrated as in Eq.([8]). The response 
function is obtained by dividing the resulting longitudinal cross section by the corresponding factors 
evaluated at the effective photon momentum, GE(g 2 yj — w 2 ) ^\—uJ 2 /q 2 j^j / {q 2 ^—uj 2 ). In the EMAr 
analysis, there is an exact cancellation of these factors. 




FIG. 5: Response function for the Is shell based on the PWIA, EMAr and full 6D calculations for 208 Pb 
at E, = 500 Mev and q = 550 MeV/c. 

The DWBA and EMAr results agree very well for the shift of the peak response relative to the 
peak of the PWIA response function. This is confirmed by fits of each of these responses to the 
EMA form of Eq. ([3]) using the parameters shown in Tabled The shift of the ls-shell peak response 
relative to the PWIA peak response is fit by 6k of -22.0 MeV for the full DWBA case (labeled 
6D-F 

exact in the table) and by -22.0 MeV for the EMAr response function. The magnitude of the 
peak EMAr response function is about 1.5% lower than that for the DWBA response function as 
shown by the A fit parameters in the table. 

In order to test one of the assumptions of the EMA or EMAr analysis, we also calculated full 
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TABLE I: EMA fit parameters for the Is shell response functions at E=500 MeV and g=550 MeV/c. 5k is 
in MeV. 





6k 


Jema 


A 


EMAr 


-22.0 


0.80 


0.985 


6-D — F exac t 


-22.0 


0.80 


1.00 


6D - F eff 


-22.0 


0.80 


1.00 



DWBA cross sections with the form factor and current-conservation factor evaluated at q e ff and 
removed from the integral over the photon momentum. However, the photon propagator was left 
within the integral. This allows a comparison of results for g^ DWBA ^> Wi th the nonlocality of the 
form factor integrated over versus otherwise identical results with the form factor evaluated at 
the effective momentum transfer and taken out of the integral. In the latter case the form factor 
simply cancels out of the result. When g( BWBA ^ j s calculated with the form factor evaluated at 
q e ff and factored out of the integral, thus canceling, the result is fit using the parameters in the 
line of Table [J labeled 6D-F e yy. When g)P WBA ) [ s calculated with the form factor integrated over 
and then divided out at q e ff, the result is fit using the parameter values in the line labeled 6D- 
F "exact in the table. Both calculations are found to produce essentially the same response functions, 
S£ WBA , in the sense that the same fitting parameters describe both equally well. Thus, there is 
not any evidence for errors associated with evaluating the form factor at q e ff and factoring it out 
of the integral. This is a nontrivial and important result because the form factor reduces the cross 
sections by about a factor 4 for q=0.55 GeV/c. There would be a significant difference if the form 
factor evaluated at the momentum transfer of the electron were divided out in Eq, (|3ip as has been 
assumed to be the correct procedure in some works. Use of the momentum transfer of the electron, 
q, versus the effective photon momentum, q e ff, leads to a difference in cross sections by a factor 
G 2 E (q 2 - u 2 )/G%(ql ff - uj 2 ) m 1.23 for 208 Pb at q = 0.55 GeV/c and u = 0.17 GeV. We find very 
clear evidence from this analysis that the form factor should be evaluated at the effective photon 
momentum transfer rather than the momentum transfer of the electron when response functions 
are extracted from data. 

We draw the following conclusions from these tests. The EMAr approximation provides a good 
approximation to the full 6D DWBA analysis. It reproduces the shift of the DWBA response 
relative to the PWIA response very well, i.e., the momentum shift 5k is the same: 
5k( EMAr '}). When the focus factors are kept within the integration over r as in the EMAr analysis, 
they do not cancel precisely with the prefactor Q^/Q^^. The r-integration provides a normalization 
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reduction of about 1.5%, i.e., A = 0.985 in fits of the EMAr results to the EMA form. In the 6D 
analysis with the nonlocality of the photon propagator also included in the integration over photon 
momentum, but everything else the same as in the EMAr calculation, there is no normalization 
correction, i.e., A = 1.00 in fits to the EMA form. We conclude that the nonlocality of the photon 
propagator produces a normalization 1.5% greater than the normalization of the EMAr response 
function: ( R^ D ^ ~ 1.015 R ^ r ' ), thus canceling the normalization reduction of the EMAr 
result. In order to include the nonlocality of the photon propagator, the normalization of the EMAr 
response for 208 Pb should be increased by the factor 1.015. It then agrees with the normalization of 
the full DWBA result because the nonlocality of the photon propagator cancels the reduction that 
arises in the EMAr result. The renormalized EMAr result is found to give excellent agreement with 
the full DWBA results for both the shift and the normalization. When the nucleon form factor 
and current-conservation factors also are kept in the integration over photon momentum, there is 
no additional change of the normalization compared with evaluating those factors at q e ff- 

D. Comparison of EMAr and EMA calculations 



The three-dimensional integral of Eq. (|17p is dominated by a stationary phase point that may 
be obtained by approximating the eikonal phase x( r ) ~ x(0) + r • Vx(0) + • • • . The effective 
momentum is then 

q e // = q + v x (o) (32) 

and the integral for the time-component of the current takes the form [using h® = cos^9 e ] 

KL MAr ' = coa^e*® J d 3 re l (^//-P)' r / l (r)//(r)^ m (r). (33) 

Generally it is found that the use of q e ff overestimates the Coulomb corrections unless Vx(0) is 
reduced by a factor in order to simulate an average value over the nucleus, i.e., 

q e // = q + /£AfAVx(0). (34) 

We refer to this stationary-point analysis with the full r-dependence of the focus factors left within 
the integral as EMAr' and use Jema = 0.8, as is consistent with fits of the EMAr result. When 
the focus factors are also approximated using 

/i(r)//(r) « (l - fEMAV(0)/ki) (l - f E MAV(0)/k f ) 
_ kj,effkf, e ff 

kikf 
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then they are cancelled by the \jQ 2 e ^ factor in the response function. That approximation leads 
to the usual EM A result, 

M nlm -COS-fJ e — WnlmK^eff ~ PJ • W 

i f 

The effective momentum approximations provide a good reproduction of the full 6D analysis for 
both the longitudinal response function, Rl, and the transverse response function, Rt, as shown in 
Figures O and [3 The EMAr' stationary-point analysis using Jem A = 0.8 produces essentially the 
same results as the EMAr that includes the integration over the variation of x( r )- The EMA result 
is also very close to the results based on EMAr and EMAr'. Thus, it is clear that the integration 
over r that is incorporated in the EMAr analysis provides results that differ only in the fine details. 
The usual EMA analysis is almost as good once one has in hand a reasonable value of Jem A to use. 
Some numerical values are given in Table |TT] in order to provide a more quantitative comparison of 
the approximations with the 6D calculation. 




co (GeV) 



FIG. 6: Longitudinal response function for the ls-shell of 208 Pb based on the PWIA (dotted line), EMA (o 
symbols), EMAr' (dash line), EMAr (solid line) and full 6D calculations (x symbols) at Ej = 500 Mev and 
q = 550 MeV/c. 
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FIG. 7: Transverse response function for the ls-shell 208 Pb based on the PWIA (dotted line), EMA (o 
symbols), EMAr' (dash line), EMAr (solid line) and full 6D calculations (x symbols) at E,; = 500 Mev and 
q = 550 MeV/c. 

TABLE II: Numerical values for R L and R T for the ls-shell of 20S Pb at E=500 MeV and g=550 MeV/c. 





R {PWIA) 


R (EMA) 


R (EMAr) 


R (EMAr') 


4 6D) 


0.10 


0.47 


0.15 


0.12 


0.15 


0.13 


0.16 


6.5 


5.4 


5.3 


5.5 


5.2 


0.20 


2.5 


3.9 


4.1 


3.9 


3.9 


to 


ltrp 




R {EMAr) 


R {EMAr') 




0.10 


0.33 


0.11 


0.09 


0.11 


0.09 


0.16 


4.6 


4.1 


4.1 


4.2 


3.9 


0.20 


1.9 


3.1 


3.2 


3.1 


3.1 



IV. CROSS SECTION CALCULATIONS 

Numerical calculations using the EMAr analysis have been performed including all shells of 
shell- model wave functions with the harmonic oscillator parameter adjusted so that the correct 
nuclear charge radius is obtained. Table HTT1 shows the parameter values used. Harmonic oscillator 



21 



TABLE III: Parameters used in calculations: (3 is the harmonic oscillator parameter; Vq and R are the 
Coulomb potential parameters. 



Nucleus 


/3(fm) 


V (GeV) 


R (fm) 


208 Pb 


3.564 


0.0256 


7.10 


56 Fe 


2.854 


0.0124 


3.97 



wave functions are used for the shell model. In coordinate space they are 

Vw(r) = NY lm (n r )r l ^(-(n - l)/2, 1 + 3/2, ryft/p) e^'^ (37) 

with normalization constants N determined by J d 3 r\ip(r)\ 2 = 1. Furthermore, Y\ m are the well 
known spherical harmonics and \F\ the confluent hyper geometric functions. 

Neutron contributions to cross sections are required in order to include the magnetic scattering. 
They are assumed to be proportional to the proton contributions and are included by using suitable 
form factors, i.e., 



N 



G 2 m > G 2 Mp + —G 2 Mn (38) 

times the proton contributions, where subscripts p and n refer to the proton and neutron, respec- 
tively. Dipole form factors 1/(1 + Q 2 /0.71GeV 2 ) 2 are used for the variation of F\ and F2 with 
Q 2 . 

Figure [8] shows cross sections for 208 Pb at 500 MeV electron energy. The momentum transfer is 
held fixed at 550 MeV/c and therefore the scattering angle varies with energy loss ui from about 70° 
to about 100°. The figure shows the plane- wave (PWIA) cross sections as light lines and the EMAr 
results for Coulomb distorted cross sections as heavy lines. Generally the Coulomb corrections shift 
the peaks to larger energy loss. In our calculations, the average binding energy of 8 MeV was used 
for all shells. Figure shows similar cross sections at 800 MeV electron energy with the momentum 
transfer held fixed at 900 MeV/c. In this case the scattering angle varies from about 70° to about 
105°. At the higher momentum transfer, the longitudinal cross section is seen to be a small fraction 
of the total cross section even without the pionic contributions. Longitudinal response functions 
for 208 Pb are shown in Figure [TUl In this figure, the PWIA response functions shown obey the 
Coulomb sum rule in the form 

J 1 + 4M{M+ui 
22 
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FIG. 8: Ratio of EMAr and PWIA longitudinal, transverse and total cross sections to Z times the Mott 
cross section for 500 MeV electron scattering from 208 Pb at q=550 MeV/c. 

where the kinematical factor in the denominator cancels the kinematical factors due to wave func- 
tions normalizations and currents. The correction is modest: the denominator factor is about 1.01 
at uj = 0.2 GeV and 1.03 at u = 0.4 GeV. 

Note that the contributions owing to correlations in the nuclear wave functions are omitted in 
our calculations. If they are small at the q values shown, the Coulomb sum rule should be satisfied 
approximately. Our wave functions are approximate and final-state interactions of the knocked-out 
nucleon have been omitted. Cross sections presented in this work are not expected to be very close 
to experimental results, however, the nuclear model used is expected to be adequate for testing the 
accuracy with which Coulombic effects can be removed. 

We have fit the EMAr longitudinal response functions to the EMA form as in Eq. ([3]) using the 
same value of q e ff in the prefactor of Eq. (|30|) as in the PWIA response function. The fits of the 
EMAr response function based on all shells yield similar values for A and a little smaller values 
for <5k compared with fits of the Is shell response function. The fit parameters are summarized 
in Table IIVI Accounting for the nonlocality of the photon propagator as in the results based on 
the 208 Pb ls-shell DWBA response, the A EMAr factors are renormalized by the factor 1.015 to 
estimate factors A DWBA for a full DWBA analysis that includes all shells. Our results support the 
use of the EMA fits of experimental data as in Eq. ([3]) using A = 1. The Jem A factors are a little 
smaller when all shells are included. That is understandable because higher shells include ones 
with wave functions that vanish at r = 0. For those shells, the distortion effects in the electron 
waves contribute at radii away from r = where the Coulomb potential is weaker. The results for 
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co (MeV) 

FIG. 9: Ratio of EMAr and PWIA longitudinal, transverse and total cross sections to Z times the Mott 
cross section for 800 MeV electron scattering from 208 Pb at q=900 MeV/c. 

TABLE IV: EMA fit parameters for the EMAr response functions at energy E in GeV and momentum 
transfer q in GeV/c. The value of 8k is in MeV. 



Nucleus 


E 


q 


5k 


Jem a 


j^EMAr 


j^DWBA 


208 Pb 


0.5 


0.55 


-21.0 


0.82 


0.98 


wl.00 


208 Pb 


0.8 


0.90 


-19.5 


0.76 


0.985 


wl.00 




0.5 


0.55 


-8.8 


0.71 


0.99 


wl.OO 


58 Fe 


0.8 


0.90 


-9.5 


0.77 


1.00 


wl.OQ 



the Is shell show that both EMAr and DWBA yield the same value of 5k. The momentum shifts 
should be equal also for response functions based on the sum over all shells. 

Figures [IT] and [12] show the longitudinal and transverse response functions Rl and Rt that 
do not have form factors divided out for Pb at q=550 MeV/c. Results are shown for the PWIA, 
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FIG. 10: EMAr longitudinal response functions divided by Z and the nucleon form factor for 208 Pb at Ej 
= 500 Mev and q = 550 MeV/c (solid line) and at E* = 800 Mev and q = 900MeV/c (dashed line). The 
corresponding PWIA response functions without Coulomb effects included are shown by the light dashed 
lines. Fits of the response functions using Eq. ([3]) are shown by the x symbols and the parameters of the 
fits are given in Table HVl 

EMA and EMAr calculations, where the EMA results are not a fit but rather are a straightforward 
calculation using Jema = 0.8 and A = 1.0. The Coulomb effects of the EMAr analysis are well 
approximated by the EMA calculation. As has been discussed, the overall magnitude of the EMA 
response is higher than the EMAr response by about 2% because the A parameter has not been 
used. 

Figure fl3l shows cross sections for 5e Fe at 500 MeV electron energy and q= 550 MeV/c and 
Figure [TH shows cross sections at 800 MeV electron energy and q = 900 MeV/c. Coulomb effects 
are somewhat smaller for the 56 Fe nucleus because the Coulomb potential is smaller. Response 
functions for the 56 i ? e target are shown in Figure [T5l Fits of the response functions to the EMA 
form of Eq. @ yield the fitting parameters shown in Table HVl The shifts are given in this case 
by Iema = 0.71 (500 MeV) and 0.77 (800 MeV). 
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FIG. 11: EMAr longitudinal response functions divided by Z for 208 Pb at Ej = 500 Mev and q = 550 
MeV/c (solid line). The corresponding PWIA response function without Coulomb effects included is shown 
by the dashed line and the EMA calculation using Jem A = 0.8 and A — 1.0 is shown by the x symbols. 

V. CONCLUSIONS 

In this work we have tested some assumptions that have been used in the analysis of experimental 
data for quasielastic scattering from nuclei. The main focus is to use a known nuclear model (in 
this case the shell model) in order to test how well the Coulomb corrections can be removed 
from DWBA cross sections using the effective- momentum approximation (EMA). The goal is to 
extract PWIA response functions from the DWBA cross sections. It is assumed that the Coulomb 
corrections are not much affected by the nuclear model used. 

At the electron beam energies considered in the work, namely 500 MeV and higher, the Coulomb 
effects in quasielastic scattering from nuclei can be described accurately using the eikonal distorted 
waves that include higher-order corrections. The eikonal analysis has simplifying features because 
one can isolate the phases that cause shifts of the electron momenta, the focusing factors and the 
spin phases that affect the L/T separation. We have used the analytical phases up to order 1/k 2 in 
the eikonal expansion that were developed in Ref. [37]]. As one check on the numerics, the eikonal 
phases were computed two ways: by direct numerical integration of the defining equations and by 
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FIG. 12: EMAr transverse response functions divided by Z for 208 Pb at Ej = 500 Mev and q = 550 MeV/c 
(solid line). The corresponding PWIA response function without Coulomb effects included is shown by the 
dashed line and the EMA calculation using Jema — 0.8 and A = 1.0 is shown by the x symbols. 

use of the analytical formulas. Both give the same results. For the cases considered in this work, 
the eikonal wave functions provide very well converged results. As a check of the three-dimensional 
integration used in the EMAr analysis, the PWIA results were computed two ways: using analytical 
Fourier transforms of the nucleon's bound state wave functions and by three-dimensional numerical 
integration. The latter calculations are the same as those for the EMAr amplitude except that the 
Coulomb effects are omitted. With suitable integration grids the results are essentially the same 
at an accuracy better than 1% near the peak of response functions and errors at larger co can be 
1% or 2% of the peak value of the response function. Generally the errors in numerical results are 
insignificant in the plots. 

Full DWBA computations are extremely time consuming. An approximation called EMAr is 
used to simplify the analysis. The EMAr analysis evaluates the full r-dependence of the eikonal 
distorted waves but approximates the hard-photon propagator and the form factor in the nucleon 
current by evaluating them at the effective momentum, q e ff- Tests of the EMAr against the full 
DWBA analysis were carried out for the response function of the Is shell of 208 Pb. Those tests 
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FIG. 13: Ratio of EMAr and PWIA longitudinal, transverse and total cross sections to Z times the Mott 
cross section for 500 MeV electron scattering from 56 Fe at q=550 MeV/c. 

showed that the EMAr produces close agreement with the DWBA. Moreover, the assumption 
that one should remove the nucleon form factor (which is integrated over in the DWBA analysis) 
by evaluating it at q e ff was found to be justified with better than 1% accuracy. This should 
be compared with large differences in cross sections when the form factor is evaluated at q, the 
momentum transfer of the electron, instead of q e ff- We find clear evidence that the form factor 
should be evaluated at the effective momentum when it is divided out of experimental cross sections 
in order to check the Coulomb sum rule. 

The analysis of Bates experimental data in Ref. [{J] uses the form factor at q rather q e ff for a 
40 Ca nucleus. Results for the Coulomb sum rule are about 0.8-0.9 compared with the expectation 
of unity. If q e fj were used in the analysis, the Bates results for the Coulomb sum rule would be 
increased by about 5%, thus making them closer to unity. The analysis of Saclay experimental 
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FIG. 14: Ratio of EMAr and PWIA longitudinal, transverse and total cross sections to Z times the Mott 
cross sections for 800 MeV electron scattering from 5G Fe at q=900 MeV/c. 



data in Refs. 
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121 ] uses form factors at q e ff- Significantly lower values for the Coulomb sum rule 
are found based on the Saclay analysis. The differences between the Bates and Saclay results are 
much larger than can be attributed to Coulomb corrections. 

We find that the spin phases in electron wave functions produce very small effects at energies 
of 500 MeV or higher. The helicity matrix elements that involve the spin phases are very close 
to those of a PWIA analysis for quasi-elastic scattering. Consequently, the Rosenbluth separation 
extracts response functions R e ^ pt and R^ pt that are accurate in the sense that they correspond 
very closely to the distorted wave matrix elements of the longitudinal and transverse parts of the 
currents. 

The effects of the distorted waves on the longitudinal response function are twofold: 1.) for 
electron scattering they shift the peak of the response functions towards larger values of the energy 
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FIG. 15: EMAr longitudinal response functions divided by Z and the nucleon form factor for 56 Fe at E, 
= 500 Mev and q = 550 MeV/c (solid line) and at E* = 800 Mev and q = 900MeV/c (dashed line). The 
corresponding PWIA response functions without Coulomb effects included are shown by the light dashed 
lines. Fits of the response functions using Eq. ([3]) are shown by the x symbols and values of the fitting 
parameters are given in Table IIVI 

loss, u), and 2.) they distort the shapes of the response functions, more so for the inner shells 
than the outer ones. However, reasonably accurate fits of the distorted response functions can be 
obtained using the EMA fitting procedure of Eq. (|3j). The momentum shift parameter 5k is found 
to be given by /em A « 0.80, for both the 2m Pb and 56 Fe nuclei, i.e., 5k « 0.8014(0), where V c (r) 
is the Coulomb potential. More precise values are given in Table HVl The normalization parameter 
A is equal to 1.00 within one or two percent. The uncertainty arises because the normalization 
for the sum over shells has been calculated based on the full DWBA for the ls-shell of 208 Pb and 
because the shape of the distorted response function differs a little from the shape of the PWIA 
response function for oj significantly away from the peak. Therefore fits to the PWIA shape cannot 
reproduce the response precisely. Note that the good agreement of /em A and A for 208 Pb and 56 Fe 
demonstrates that the Coulomb corrections do not depend significantly on the nuclear model. Note 
also that the analysis of experimental data using a fit as in Eq. ([3|) tends to give more accurate 
results for pS^ WIA ^ at the peak of the response because that is controlled by 5k and less accurate 
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results away from the peak because of the distortion of the shape. 

Estimates of longitudinal, transverse and total cross sections have been calculated using shell 
model wave functions for 208 Pb and 56 Fe at q = 0.55 GeV/c and q = 0.8 GeV/c. These kinematical 
conditions match the ones used in a recent experiment at Jlab. Because final-state interactions, 
correlations and pion production have been omitted, the calculated cross sections may differ signif- 
icantly from experimental cross sections. Nevertheless the Coulomb corrections should be reliable 
at the level of a few percent. 

Coulomb corrections are notoriously difficult to calculate and our calculations refute claims 

□ 

that may be found in the literature. For example, Ref. [461 ] claims that the EMA procedure is not 
accurate for the longitudinal response at 485 electron energy and 60° scattering angle for a 208 Pb 
target. The basis for the claim is that significant differences are found between EMA results and 
results based on an ad-hoc DWBA analysis that has been used extensively. We find that the EMA 
with appropriate parameters can describe the ls-shell DWBA or all-shells EMAr results very well 
at essentially the same kinematics. We wish to emphasize that all of our numerics are under good 
control and various consistency checks have been made that give confidence in the results reported 
herein. 
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